Parity violation in atomic ytterbium: experimental sensitivity and systematics 



K. Tsigutkin/'B D. Dounas-Frazer, 1 A. Family, 1 J. E. Stalnaker, 1 ^ V. V. Yashchuk, 2 and D. Budker 1 ' 3 

'Department of Physics, University of California at Berkeley, 
Berkeley, CA 94720-7300 
2 Advanced Light Source Division, Lawrence Berkeley National Laboratory, Berkeley CA 94720 
3 Nuclear Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720 

(Dated: January 4, 2010) 

We present a detailed description of the observation of parity violation in the 6s 2 1 So — V 5d6s 3 Di 
408- nm forbidden transition of ytterbium, a brief report of which appeared earlier. Linearly polarized 
408-nm light interacts with Yb atoms in crossed E- and B-fields. The probability of the 408- 
nm transition contains a parity violating term, proportional to {£ ■ B)[(E X £ ) • B], arising from 
interference between the parity violating amplitude and the Stark amplitude due to the E-field (£ 
is the electric field of the light). The transition probability is detected by measuring the population 
of the 3 Po state, to which 65% of the atoms excited to the 3 Di state spontaneously decay. The 
population of the 3 Po state is determined by resonantly exciting the atoms with 649-nm light to 
the 6s7s 3 Si state and collecting the fluorescence resulting from its decay. Systematic corrections 
due to E-field and B-field imperfections are determined in auxiliary experiments. The statistical 
uncertainty is dominated by parasitic frequency excursions of the 408-nm excitation light due to 
imperfect stabilization of the optical reference with respect to the atomic resonance. The present 
uncertainties are 9% statistical and 8% systematic. Methods of improving the accuracy for the 
future experiments are discussed. 

PACS numbers: 11.30.Er, 32.90.+a 



I. INTRODUCTION 

In an earlier paper we reported on observa- 
tion of the atomic parity violation (PV) effect in the 
6s 2 1 So — > 5d6s 3 Di 408-nm forbidden transition of 
174 Yb. We measured the PV induced transition matrix 
element to be 8.7 ± 1.4 x 10~ 10 cao, which confirms the 
theoretically anticipated PV enhancement in Yb |2| and 
constitutes the largest atomic parity violation effect ob- 
served so far. However, the measurement accuracy is 
not yet sufficient for the observation of the isotopic and 
hyperfine differences in the PV amplitude, the study of 
which is the main goal of the present experiments. Here 
we describe the impact of the apparatus imperfections 
and systematic effects on the accuracy of the measure- 
ments and discuss ways of improving it. 

During the initial stage of the experiment, an effort 
was invested into measuring various spectroscopic prop- 
erties of the Yb system of direct relevance to the PV 
measurement, including determination of radiative life- 
times, measurement of the Stark- induced amplitudes, hy- 
perfine structure, isotope shifts, and dc-Stark shifts of 
the 6s 2 1 So — > 5d6s 3 Di transition Q. In addition, the 
6s 2 % -> 3 D 2 transition at 404 nm has been observed, 
and the electric quadrupole transition amplitude and ten- 
sor transition polarizability have been measured [lj . The 
forbidden magnetic-dipole (Ml) amplitude of the 408 nm 
transition was measured to be 1.33 x 10 using the 
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Ml-(Stark-induced)El interference technique [5j. The 
ytterbium atomic system, where transition amplitudes 
and interferences are well understood, has proven use- 
ful for gaining insight into the Jones-dichroism effects 
that had been studied in condensed-matter systems at 
extreme conditions and whose origin had been a matter 
of debate (see Ref. |(| and references therein). 

An experimental and theoretical study of the dynamic 
(ac) Stark effect on the 6s 2 1 So — > 5d6s 3 Di forbid- 
den transition was also undertaken Q ■ A model was de- 
veloped to calculate spectral line shapes resulting from 
resonant excitation of atoms in an intense standing light 
wave in the presence of off-resonant ac-Stark shifts. A bi- 
product of this work was an independent determination 
(from the saturation behavior of the 408-nm transition) 
of the Stark transition polarizability, which was found to 
be in agreement with the earlier measurement 

The present Yb APV experiment involves a measure- 
ment using an atomic beam. An alternative approach 
would involve working with a heat-pipe-like vapor cell. 
Various aspects of such an experiment were investigated, 
including measurements of collisional perturbations of 
relevant Yb states [1], nonlinear optical processes in a 
dense Yb vapor with pulsed UV-laser excitation [§], and 
an altogether different scheme for measuring APV via 
optical rotation on a transition between excited states 
10]. 

The present paper addresses the issues of sensitivity 
and systematics in the Yb APV experiment. In Sec- 
tions |TT] and Mil the experimental technique and its ap- 
plication in the present experiment are discussed. In 
Section IIVI a method of analyzing the impact of various 
apparatus imperfections is described based on theoret- 
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ical modeling of signals recorded by the detection sys- 
tem in the presence of imperfections. In Section [V] a de- 
tailed description of the experimental apparatus is given, 
along with a discussion of the origins of the imperfections, 
which is followed by an account of the measurements of 
the imperfections in Section fVTl In Sections I VIII and IVIIII 
we discuss measurements and analysis of the PV ampli- 
tude and systematic effects, and ways of improving the 
accuracy of the PV measurements to better than 1% in 
order to measure the difference in the APV effects be- 
tween different isotopes and hyperfmc components. 



II. EXPERIMENTAL TECHNIQUE FOR THE 
APV MEASUREMENT 

As discussed in Ref. [l| , the idea of the experiment is 
to excite the forbidden 408-nm transition (Fig. [TJ with 
resonant laser light in the presence of a quasi-static elec- 
tric field. The PV amplitude of this transition arises 
due to PV mixing of the 5d6s 3 Di and 6s6p 1 Pi states. 
The purpose of the electric field is to provide a refer- 
ence transition amplitude, which is due to Stark mixing 
of the same states interfering with the PV amplitude. 
In such an interference method (ill Il2| , one is measur- 
ing the part of the transition probability that is linear in 
both the reference Stark-induced amplitude and the PV 
amplitude. In addition to enhancing the PV dependent 
signal, the Stark-PV interference technique provides for 
all-important reversals that separate the PV effects from 
the systematics. 




FIG. 1: (color online) Low-lying energy eigenstates of Yb and 
transitions relevant to the APV experiment. 

Even though the APV effect in Yb is relatively large, 
and the Ml transition is strongly suppressed, the Ml 
transition amplitude is still three orders of magnitude 
larger than the PV amplitude. As a result, the geometry 
of the experiment was designed to suppress spurious Ml- 



Stark interference. In addition, this effect is minimized 
by the use of a power-build-up cavity to generate a stand- 
ing light wave. Since a standing wave has no net direction 
of propagation any transition rate which is linear in the 
Ml amplitude, will cancel out (see below). 

The advantages of the present experimental configura- 
tion can be demonstrated by considering Yb atoms in the 
presence of static electric, E, and magnetic, B, fields in- 
teracting with a standing monochromatic wave produced 
by two counter-propagating coherent waves in an optical 
cavity. The electric field in the standing wave, £, is a sum 
of the fields of the two waves. For resonant atoms, the 
transition rate from the ground state 1 Sq to the excited 
state 3 Di is (see e.g. pj], Eq. (3.127)) 

Rm = j^\A m \ 2 , (I) 

where T is the natural linewidth of the transition, Am 
is the transition amplitude, and M = 0, ±1 is the mag- 
netic quantum number of the excited state. Here and 
in the rest of this section it is assumed that the individ- 
ual magnetic sublevels of the 3 Di state are resolved. For 
convenience, we set h = 1 henceforth and measure the 
transition rate in units of T. 

The transition amplitude Am is the sum of the electric- 
(El) and magnetic- (Ml) dipole transition amplitudes: 

A M = A M +A m . (2) 

The El amplitude has two contributions corresponding 
to the Stark- and PV- mixing of the 3 Di and 1 Pi states. 
That is, 

.(El) _ . (Stark) ,(APV) 
M — M M 

= ip{-i) M (e x s)_ M + K{-i) M e- M , (3) 

where /? is the vector transition polarizability, C is re ~ 
lated to the reduced matrix clement of the Hamilto- 
nian describing the weak interaction, and £o,±i are t ne 
spherical components of the vector £. Although Stark- 
induced transitions are generally characterized by scalar, 
vector, and tensor polarizabilities 0, [ll|, for the case of 
a J = — » 1 transition, only the vector polarizability 
contributes. Equation ([3]) is derived in Appendix [XJ 

Similarly, the Ml transition amplitude has two com- 
ponents: one for each of the two counter-propagating 
laser beams. Let £ + = £ + £ and £ = £_ £ denote the 
electric fields of the beams traveling in the k and — k di- 
rections, respectively. Then £ = £ + + £ and the Ml 
amplitude is given by 

A^ = M(-1) M (k x £ + )_ M + M(-l) M (-k x £_)_> 
= M(-l) M (Skx£)_ M , (4) 

where M. is the reduced matrix element of the Ml tran- 
sition and k is a unit vector in the direction of the 
wavevector. Here we have introduced <5k = 5k k with 
6k = (£ + —£_)/£. For a perfect standing wave, £+ = £- 
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and hence 5k = and the Ml transition is completely 
suppressed. In practice, £_ = £+ — 8£ due to the small 
but nonzero transmission of the back mirror in the cavity. 
Since \5£\ < £, \Sk\ m \S£/£\ < 1. Thus, although the 
Ml transition amplitude is not strictly zero, it is highly 
suppressed. 

Without loss of generality, the quantities f3, £, and M. 
are assumed to be real. In general, the rate Rm given 
by Eq. ([1]) includes terms proportional to (3M. (Stark-Mi 
interference) and f3( (Stark-PV interference). 

A careful choice of field geometry allows for additional 
suppression of undesirable Stark-Mi interference. From 
Eq. ©, it is evident that the Stark-PV interference is 
proportional to the rotational invariant 

( £ ■B)[(Exf).B], (5) 



which is P-odd and T-even. In the present experimental 
apparatus the electric field, E, is applied orthogonally to 
the magnetic field, B, and collinearly with the axis of the 
linearly-polarized standing light wave, as shown in Fig. [5] 




FIG. 2: (color online) Orientation of fields for PV-Stark in- 
terference experiment and schematic of the present APV ap- 
paratus. Not shown is the vacuum chamber containing all the 
depicted elements, except the photomultiplier (PMT) and the 
photodiode (PD). PBC-power buildup cavity. Light is ap- 
plied collinearly with x. 

This geometry is such that the Ml and Stark-induced 
amplitudes are out of phase. Thus, they do not interfere 
and therefore do not produce spurious PV-mimicking ef- 
fects (see Section |IV}. 



III. PV SIGNATURE: IDEAL CASE 

In the ideal case where we neglect the apparatus im- 
perfections, the static magnetic and electric fields are 
B = Bz and E = Bx, respectively, and the light stand- 
ing wave has an electric field 

£ = £(sin0y + cos0z). (6) 



With this field orientation (see Fig. [2]), Eqs. ((T|) through 
© yield 

Ro = 4£ 2 [p 2 E 2 sin 2 9 + 2(/3E sin cos 0) , (7) 
R±! = 2£ 2 (/3 2 E 2 cos 2 9~2(fiE sin cos 6) , (8) 

where terms of order £ 2 and higher are neglected, and 
5k = is assumed. 

In order to isolate the Stark-PV interference term from 
the dominant Stark-induced transition rate, we modulate 
the electric field: E — E dc + Eg cos(wf), where Eq is the 
modulation amplitude, w is the modulation frequency, 
and E dc provides a DC bias. Then Eqs. J7]) and © 
become 

Rm = Km + K M ] cos{ujt) + TZ M ] cos(2wi) , (9) 

In] 

where TZ M is the amplitude of the nth harmonic of the 
transition rate Rm- The dominant Stark-induced con- 
tribution, which oscillates at twice the modulation fre- 
quency, is 

TZ l2] = 2fi 2 E 2 £ 2 sin 2 9, (10) 
TZ l l\ = (3 2 E 2 £ 2 cos 2 0. (11) 

On the other hand, the amplitude IZ^j contains the 
Stark-PV interference term: 

= 8£ 2 {j3 2 EqE Ac sin 2 + CPEq sin cos 9 \ , (12) 

K l l\ = A£ 2 {p 2 EqE Ac cos 2 - CPEq sinfl cos 6>) . (13) 

The term IZ^j is a constant "background" : 

K 0] = 4£ 2 (p 2 {El + Ej c ) sin 2 + AQf3E dc sin cos 9 \ , 
U l l\ = 2£ 2 (p 2 {El + El c ) cos 2 - 4C/3£ dc sin6»cos6l) . 

For an arbitrary polarization angle 9, all three Zee- 
man components of the transition, as shown in Fig. [3^, 
are present while scanning over the spectral profile of the 
transition. The first-harmonic signal due to Stark-PV in- 
terference has a characteristic signature: the sign of the 
oscillating terms for the two extreme components of the 
transition is opposite to that of the central component. 
The second-harmonic signal provides a reference for the 
lineshape since it is free from interference effects linear 
in E (Fig. [2b) ■ With a non-zero DC component present 
in the applied electric field, a signature identical to that 
in the second harmonic will also appear in the first har- 
monic, Fig. [3fc. The latter can be used to increase the 
first-harmonic signal above the noise, which makes the 
profile analysis more reliable. 

To obtain the PV term from the measured first- 
and second-harmonic transition rates, we first normalize 
the first-harmonic signals TZ^J by their second-harmonic 
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FIG. 3: Discrimination of the PV effect by E-field modulation 
under static magnetic field. The Zeeman pattern is shown for 
the polarization angle 6 = 7r/4. 
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counterparts 7£U and combine the results in the follow- 
ing way: 



K 





K ll] 


n [1] 




K l2] 





= T 



16C 

Wo 



(14) 



Here we take 6 = ±7r/4, which are optimal polarization 
angles for the PV measurements (see next section). This 
method has the advantage that K, is independent of E& c , 
so that the E-field bias may be chosen based on technical 
requirements of the experimental apparatus. 



IV. PV SIGNATURE: IMPACT OF 
APPARATUS IMPERFECTIONS 

While the current Yb-APV apparatus has been de- 
signed to minimize systematic effects, the PV mimicking 
systematics may be a result of a combination of multiple 
apparatus imperfections. In order to understand the im- 
portance of these effects, the electric and magnetic field 
misalignments and stray fields were included in a theoret- 
ical model of the transition rates as well as the excitation 
light's deviations from linear polarization. In addition, 
we relax the assumption that 5k = and include the 
effects of the residual Ml transition. 

The quantization axis is defined along z, and following 
the ideal case model, the axis of the standing light wave 
is collinear with x. We added a small ellipticity to the 
light polarization by taking 



£ = £(y sm.6 + ze* 



(15) 



where <f> is a small phase. For \(f)\ <C 1, the ellipticity 
of the light is 2<fr sin(2#) . The electric field imperfections 
are included by taking 



E = E + E', 



where 



E = (E x + e y y + e z z) cos(ojt) 
E' = E dc x + e y y + e z z, 

are the AC and DC components of the electric field. It is 
assumed that the y- and z-components of the AC field are 
in phase with the leading oscillating E-field. The impact 



of the out-of-phase AC components was analyzed within 
a complete model of the systematics and found to be 
negligible. The AC components are due to misalignments 
of the applied E-field with respect to the light wave axis 
as well as to the quantization axis z. The DC components 
arise due to a misalignment of the DC-bias field and also 
due to stray electric fields in the interaction region. 

The magnetic-field imperfections are defined within 
the same frame of reference by taking analogously 



where 



B = B + B', 



B = b x x + b y y 
B' = b' x x + b'y 



Bz 
b'z, 



where B and B' are reversing and stray non-reversing 
magnetic fields, respectively. 

Equations |T|) through ^ apply when the quantiza- 
tion axis is along the magnetic field, thus a rotation D is 
applied to each of the vectors E, B, £, and k such that 
DB oc z. That is, we take 

B DB, E -> HE, £ -> B£, and k -> Dk, (16) 

where 



(17) 



Here the matrix @(a., n) represents a rotation about an 
axis n through angle a. The angles a x and a y are given 
by 



x x , y = (B-b' z )(b' y!X +b y , x )/B 2 



(18) 



Thus, the electric field E and the polarization vector 
£ acquire additional components after the rotation (be- 
sides, for example, e y and e y ). 

Due to the imperfections, the normalized-rate modu- 
lation amplitudes now include additional terms besides 
the Stark- and the PV effects: 



n 
n 



[i] 

M _ _ (Stark) 
[2] - TM - r M 
M 



JAPV) , (Ml) 
M ' ' M 



„(</>) 
M > 



(19) 



( Stark) 

where r M is the Stark contribution due to the DC- 



bias and the field imperfections, 

„(M1) 



„(APV) 
M 



is the PV-Stark- 



interference term, r M is the Ml-Stark-intcrference con- 
tribution, and is a contribution due to the distorted 
linear polarization of the light (which is a Stark contri- 
bution, but we explicitly separate the contribution linear 
in </>). Expressions for the lowest-order terms are sum- 
marized in the Table [I] 

The normalized-amplitude combination ()14D has been 
chosen to determine the PV asymmetry. Since the Ml 
and ellipticity terms have opposite signs for M = ±1, 
their contributions to K, cancel, while the contributions 



from r 



(APV) 
M 



add. 
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TABLE I: Lowest-order terms contributing to the normalized 
transition-rate modulation amplitudes rjf. 
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The Stark-contribution, rj^ , has several terms that 
are produced due to different imperfections and impacts 
all three Zeeman components, M = 0,±1. In order to 
determine which terms could potentially mimic the PV 
asymmetry in )C, we discriminate the PV contribution 
to JC with respect to the B-field reversal and flip of the 
polarization angle, 9. Switching to a different Zeeman 
component of the transition is also a reversal, which is in- 
corporated in the expression for the asymmetry, /C. Anal- 
ysis of the noise affecting the accuracy of PV-asymmetry 
measurements demonstrate that the highest signal-to- 
noise ratio is achieved when 9 = ±7r/4, and therefore, 
the polarization flip is a change of the polarization an- 
gle by 7r/2. Thus, the normalized-amplitudc combina- 
tion (|14[) must be determined for four different combina- 
tions of the B-field directions and light-polarization an- 
gles: £(+.B,+7r/4), /C(-£,+tt/4), K(+B, -7I-/4), and 
]C(—B, — 7r/4), so that terms having different symmetries 
with respect to the reversals can be isolated: 
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The result of this procedure is summarized in Table [TT] 

TABLE II: List of the lowest-order terms contributing to the 
asymmetry K, for \8\ — 7r/4 sorted with respect to their re- 
sponse to the reversals. IC4 corresponding to a rather long list 
of terms that are invariant with respect to all reversals, is not 
shown in the table. 
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The PV asymmetry contributing to K,\ is B-field even, 
6>-flip odd. It competes with the second-order terms that 
are a combination of the E-ficld and B-ficld alignment im- 
perfections and stray fields. Using the theoretical value 
of ( ~ 10 -9 ea EH combined with the measured 
|/3| = 2.24t°;^ x 10" 8 ea /(V/cm) 0,0, the expected 



PV asymmetry, 16(//3E, is - 4 ■ for 9 = ir/4 and 

Eq = 2 kV/cm. For a typical value of misalignments 
and "parasitic" fields, e VtZ /E and b x /B (on the order of 
10~ 3 in the present apparatus), the contribution of the 
"parasitic" terms may be up to a few percent of the total 
value of K-i . Ways of measuring the contribution of these 
imperfections are discussed in the following sections. 



V. EXPERIMENTAL APPARATUS 

The forbidden 408-nm transition is excited by resonant 
laser light coupled into the power-buildup cavity in the 
presence of the magnetic and electric fields. The transi- 
tion rates are detected by measuring the population of 
the 3 Po state, where 65% of the atoms excited to the 
3 Di state decay spontaneously (Fig. [I}. This is done 
by resonantly exciting the atoms with 649-nm light to 
the 6s7s 3 Si state downstream from the main interac- 
tion region, and by collecting the fluorescence resulting 
from the decay of this state to the 3 Pi and 3 P2 states 
and subsequently, from the decay of the 3 Pi state to the 
ground state 1 So (556 nm transition). As long as the 408- 
nm transition is not saturated, the fluorescence intensity 
measured in the probe region is proportional to the rate 
of that transition. 

A schematic of the Yb-APV apparatus is shown in Fig. 
[5] A beam of Yb atoms is produced (inside of a vacuum 
chamber with a residual pressure of ~ 3 x 10~ 6 Torr) with 
an effusive source: a stainless-steel oven loaded with Yb 
metal, operating at 500 — 600° C. The oven is outfitted 
with a multi-slit nozzle, and there is an external vane 
collimator reducing the spread of the atomic beam in the 
horizontal direction. The resulting Doppler width of the 
408-nm transition is « 12 MHz 0. 




FIG. 4: (color online) Schematic of the power buildup cavity. 

Downstream from the collimator, the atoms enter the 
main interaction region where the Stark- and PV-induced 
transitions take place. Up to 80 mW of light at the tran- 
sition wavelength of 408.345 nm in vacuum is produced 
by frequency doubling the output of a TkSapphire laser 
(Coherent 899) using the Wavetrain cw ring-resonator 
doubler. After shaping and linearly polarizing the laser 
beam, « 10 mW of the 408-nm radiation is coupled into a 
power buildup cavity (PBC) inside the vacuum chamber. 
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FIG. 5: (color online) Schematic of the optical setup. Light at 408-nm is produced by frequency doubling the output of a 
Ti:Sapphire laser (Coherent 899) using the Wavetrain cw ring resonator doubler. The laser is locked to the PBC using the 
FM-sideband technique. The PBC is locked to a confocal Fabry-Perot etalon. This scannable etalon provides the master 
frequency. The 649- nm excitation light is derived from a single- frequency diode laser (New Focus Vortex). The diode laser is 
locked to a frequency-stabilized He-Ne laser using another scanning Fabry-Perot etalon. 



The cavity was designed to operate as an asymmet- 
ric cavity with flat input mirror and curved back mirror 
with a 25-cm radius of curvature and 22-cm separation 
between the mirrors. The atomic beam intersects the 
cavity mode in the middle of the cavity, where the 1/e 2 
radius of the mode in intensity is 172 fim. The asymmet- 
ric configuration has the advantage of larger mode radius 
at the interaction position compared to a symmetric cav- 
ity. A larger mode allows us to reduce the ac-Stark shifts, 
consequently reducing the width of the 408 nm transition. 
Alternatively, the cavity can be modified to operate in the 
symmetric confocal condition where multiple transverse 
modes can be excited, thereby increasing the effective 
"mode" size. However, we were unable to obtain high 
power and stable lock for the confocal configuration. 

The cavity mirrors were purchased from Research Elec- 
tro Optics, Inc. For the flat input mirror the transmission 
is 350 ppm with the absorption and scattering losses of 
150 ppm total at 408 nm. The curved back mirror is de- 
signed to have a lower transmission of 50 ppm in order 
to additionally suppress the net light wave vector and, 
therefore, the MI transition amplitude. The absorption 
and scattering losses in the curved mirror are 120 ppm. 
The finesse and the circulating power of the PBC are 
up to T — 9000 and P = 8 W. These parameters were 



routinely monitored during the PV measurements. De- 
tails of the characterization of the PBC are addressed in 
Appendix [Bl 

We found that the use of the 408-nm-PBC in vacuum 
is accompanied by substantial degradation of the mirrors. 
Typically after 6 hours of operation, the finesse drops by 
a factor of two. This is a limiting factor for the duration 
of the measurement run. The degradation of the finesse 
is due to the increased absorption and scattering losses. 
This effect is reversible: the mirror parameters can be 
restored by operating the PBC for several minutes in 
air, which makes performing a number of runs possible 
without replacing the mirrors. However, it takes several 
hours with the present apparatus to reach the desired 
vacuum after exposing the PBC to air. Presently, this 
effect is under investigation aiming for longer-duration 
experiments and shorter breaks in between. 

A schematic of the PBC setup is presented in Fig. |U 
The mirrors are mounted on precision optical mounts 
(Lees mounts) with micrometer adjustments for the hor- 
izontal and vertical angles and the pivot point of the 
mirror face. The mirror mounts are attached to an Invar 
rod supported by adjustable table resting on lead blocks. 
The input mirror is mounted on a piezo-ceramic trans- 
ducer allowing cavity scanning. 
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FIG. 6: (color online) Schematic of the E-field electrodes as- 
sembly, and a result of the E-field modeling showing an X-Z 
slice of the amplitude of the E-field z-component, e z , in a 
midplane (Y=0) of the assembly normalized by the total E- 
field amplitude, E. The voltage is applied to electrode 1, and 
electrode 2 and the correction electrodes are grounded. 



The laser is locked to the PBC using the FM-sideband 
technique pl|. In order to remove frequency excursions 
of the PBC in the acoustic frequency range, the cavity 
is locked to a more stable confocal Fabry-Perot etalon, 
once again using the FM-sideband technique. This stable 
scannable cavity provides the master frequency, with the 
power-build-up cavity serving as the secondary master for 
the laser. A schematic of the optical system is presented 
in Fig. El 

The magnetic field is generated by a pair of rectangu- 
lar coils designed to produce a magnetic field up to 100 G 
with a 1% non- uniformity over the volume with the di- 
mensions of 1 x 1 x 1 cm 3 in the interaction region. Addi- 
tional coils placed outside of the vacuum chamber com- 
pensate for the external magnetic fields down to 10 mG 
at the interaction region. The residual B-ficld of this 
magnitude does not have an impact on the PV measure- 
ments since its contribution is discriminated using the 
field reversals. 

The electric field is generated with two wire- frame elec- 
trodes separated by 2.1 cm (see Fig. [5]). The copper elec- 
trode frames support arrays of 0.2-mm diameter gold- 



plated wires. This design allows us to reduce the stray 
charges accumulated on the electrodes by minimizing the 
surface area facing the atomic beam, thereby minimizing 
stray electric fields. AC voltage of up to 10 kV at a 
frequency of 76.2 Hz is supplied to the electrodes via a 
high-voltage amplifier. An additional DC bias voltage of 
up to 100 V can be added. 

The result of the electric field non-uniformity calcula- 
tions is shown in Fig. These calculations demonstrate 
that errors in the centering of the light beam with respect 
to the E-field plates may induce substantial parasitic 
components as large as, for example, \e z \ ~ 5 x 10~ 3 E , 
producing parasitic effects comparable to the PV asym- 
metry. In order to measure and/or compensate the 
impact of the parasitic fields, additional electrodes de- 
signed to simulate stray E-field components were added 
to the interaction region. By applying high-voltage to 
these electrodes ("correction electrodes" in Fig. [6j, the 
parasitic-field components may be exaggerated and ac- 
curately measured as described in the following sections. 

Light at 556 nm emitted from the interaction region 
is collected with a light guide and detected with a pho- 
tomultiplier tube. This signal is used for initial selec- 
tion of the atomic resonance and for monitoring pur- 
poses. Fluorescent light from the probe region is col- 
lected onto a light guide using two optically polished 
curved aluminum reflectors and registered with a cooled 
photodctcctor (PD). The PD consists of a large-area 
(lxl cm 2 ) Hamamatsu photodiode connected to a 1-Gil 
transimpedance pre-amplifier, both contained in a cooled 
housing (temperatures down to — 15°C). The pre-amp's 
bandwidth is 1 kHz and the output noise is ~ 1 mV 
(rms). The 649- nm excitation light is derived from a 
single-frequency diode laser (New Focus Vortex) produc- 
ing ps 1.2 mW of cw output, high enough to saturate 
the 6s6p 3 Po — > 6s7s 3 Si transition. Due to the sat- 
uration of this transition, ~3 fluorescence photons per 
atom exited to the 3 Po state are emitted at the probe 
region. The natural width of the 649-nm transition is 
70 MHz, thus, its profile covers all transverse velocity 
groups (v x ) in the atomic beam (w 8 MHz Doppler width 
at 649 nm). A drift of the laser frequency (~ 100 MHz 
per minute) is eliminated by locking the diode laser to a 
frequency-stabilized He-Ne laser using a scanning Fabry- 
Perot etalon with the scanning rate of 25 Hz. The spec- 
tral distance between the etalon transmission peaks from 
the two lasers is measured during each scan and main- 
tained constant within an accuracy of ±3 MHz, good 
enough to eliminate any degradation of the probe-region 
signal. 

The signals from the PMT and PD are fed into lock- 
in amplifiers for frequency discrimination and averaging. 
A typical time of a single spectral-profile acquisition is 
20 s. The signals at the first and the second harmonic of 
the electric-field modulation frequency are registered si- 
multaneously, which reduces the influence of slow drifts, 
such as instabilities of the atomic-beam flux. The mod- 
ulation frequency is limited by several factors. Thermal 
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distribution of atomic velocities in the beam causes a 
spread (of « 2 ms) in the time of flight between the in- 
teraction region and the probe region. This, along with 
the finite bandwidth of the PD, leads to a reduction of 
the signal-modulation contrast (see below). The choice 
of the modulation frequency of 76.2 Hz is a tradeoff be- 
tween this contrast degradation and the frequent E-field 
reversal. 



VI. RESULTS AND ANALYSIS 

In Fig. [7] a profile of the B-field-split 408- nm spectral 
line of the 174 Yb is shown. The 649 -nm-light-induced flu- 
orescence was recorded during a single profile scan. Sta- 
tistical error bars determined directly from the spread of 
data are smaller than the points in the figure. The pecu- 
liar asymmetric line shape of the Zccman components is 
a result of the dynamic Stark effect Q ■ During a typical 
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FIG. 7: (color online) A profile of the B-field-split 408-nm 
spectral line of 174 Yb recorded at 1st- and 2nd-harmonic of 
the modulation. Also shown is a simulated PV contribution in 
the first-harmonic signal. E=5 kV/cm; DC offset=40 V/cm; 
6 — 7r/4; the effective integration time is 200 ms per point. 

experimental run 100 profiles are recorded for each com- 
bination of the magnetic field and the polarization angle 
(400 profile scans in total). In order to compute the nor- 
malized amplitude, r q , of a selected Zeeman component, 
the actual first-harmonic signal near the Zeeman peak 
is divided by the respective second-harmonic signal and 
then averaged over a number of the data points 1 . Then, 
the combination JC of Eq. (TTJ| is computed for each pro- 
file scan followed by averaging the result over all the scans 
at a given B-9 configuration. This procedure is repeated 
for all four reversals, and all B-9 symmetrical contribu- 
tions, /Ci_4, are determined. In the present experiment, 
the values of /C2,3,4-terms are found to be consistent with 



1 In the normalized rate calculations only data points having in- 
tensity higher that 1/3 of the respective Zeeman peak are used 
to avoid excessive noise from spectral regions with low signal 
intensity. 




zero within the statistical uncertainty, which is the same 
as that of the PV-asymmetry (see below). 

As can be seen from Table HH terms in K,\ associated 
with the fields imperfection are of crucial importance: 

16 

In order to measure the contribution of these terms, arti- 
ficially exaggerated E-field imperfections both static and 
oscillating, e ex , ej x , e y x and e ex , are imposed by use of 
the "correction electrodes" (see Fig. [6]), and two sets of 
the experiments were performed. In the first one, a DC- 
voltage was applied to the correction electrodes, and the 
measurements were done reversing e y x and e e z . These 

experiments yield values of e y and e z + 2Eob x / B. In the 
second set, an AC- voltage modulated synchronously with 
the main E-field is applied to the correction electrodes. 
In order to reverse the sign of the parasitic terms a tt- 
phase-shift of e y x and e ex with respect to the modulation 
signal is employed by switching the wiring of the correc- 
tion electrodes. Thus, values of the DC-imperfections, e y 
and e z , are determined. The magnitudes of the applied 
electric fields and their distributions are calculated using 
a 3D-numerical-model of the interaction region. The re- 
sults of the experiments are presented in Table UTTJ The 



TABLE III: Results of measurements of the electric field 
imperfections using artificially exaggerated AC- and DC- 
components, &y* z and e^ z . These fields were generated by 
use of the correction electrodes, Fig. [6] Eq = 2000(2) V/cm. 


DC-Set 


AC-Set 


Exaggerated imperfections (V/cm) 


eT = -140(2) 
eT = 20(2) 


eT = -120(2) 
ef = 30(2) 


Measurements (mV/cm) 


-ex 


Sb S =i6(io) 

(2S |+e,)C= 442(10) 
o 2E 




Parasitic fields (V/cm) 
ey = 3.2(2) 


e y = 0.5(0.6) 


(2£ % +e z ) = -12.6(0.3) 

JD 


e z = -1.3(0.2) 



net contribution of these imperfections to K,\ in the ab- 
sence of the exaggerated fields is found to be 2 : 



-y 




■v 



2E L> j ' 2E 
-2.6(1.6) stat .(l-5) sys t. mV/cm. 



(21) 



Compare with the PV asymmetry parameter f//3 40 mV/cm. 
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The systematic uncertainty comes from a sensitivity of 
the numerical model of the E-field, which is used for cal- 
culating the exaggerated fields in the interaction region, 
to an imperfect approximation of the electrode-system 
geometry. These experiments suggest that this field's 
imperfection cannot mimic the PV-effect entirely, never- 
theless, it appears to be a major source of systematic un- 
certainty impacting the accuracy of the PV-asymmetry 
measurements. The most prominent contribution is given 
by a combined effect of the parasitic components of the 
electric field and the non-zero projection of the lead- 
ing magnetic field on the direction of the electric field: 
e y (e z /2E + b x /B). The PV-asymmetry parameter, (/(3 
is obtained from the measured value of K\ by compensat- 
ing for the influence of these magnetic- and electric-field 
imperfections, Eq. ([2"T]) . 

There is another effect that cannot, by itself, mimic 
the PV-asymmetry, but needs to be taken into account 
for proper calibration. This effect is related to the E-ficld 
modulation implemented in the present experiment. The 
atoms are excited to the metastable state, 6s6p 3 Po, by 
the light beam in the interaction region and then travel 
^20 cm until they are detected downstream in the probe 
region. Due to the spread in the time-of-flight between 
the interaction and probe regions, the phase mixing leads 
to a reduction of the signal modulation contrast at the 
probe region, and it depends on the signal-modulation 
frequency. Since the signal comprises two time-scales of 
interest, first- and second-harmonic of the E-field mod- 
ulation, the contrast reduction is different for the two. 
Therefore, the ratio of the signal modulation amplitudes, 
TMi on which we base the PV-asymmetry observation, 
appear altered in the probe region compared to what it 
would be at the interaction region. The amplitude com- 
bination, /C, and, therefore, the PV-parameter, are 
similarly affected. In our data analysis, a correction co- 
efficient, Co, is introduced, which has been calculated 
theoretically: 
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FIG. 8: (color online) The PV interference parameter £//3. 
Mean value: 39(4) atat . (3) syst . mV/cm, |C| = 8.7 ± 1.4 x 
10~ 10 ea . 



model of the transition rates employing the field geome- 
try shown in Fig. [51 The direction and, thus, the signs of 
the electric and magnetic fields as well as the polarization 
angle 9 were calibrated prior to the PV measurements. 
Special care was taken of detecting parasitic phase shifts 
in the lock-in amplifier. A signal from an arbitrary func- 
tion direct digital synthesis (DDS) generator simulating 
the output of the probe region photodetector was fed 
into the amplifier. The signal is comprised of a sum of 
two sinusoidal waveforms, one frequency doubled, attenu- 
ated, and phase shifted with respect to the other. Results 
of the signal parameters measurement from the lock-in, 
such as the first-to-second harmonic amplitude ratio, rel- 
ative phase shift and its sign, are compared to those used 
in the DDS generator to simulate the signal. The differ- 
ence in the measured and generated amplitude ratio is 
found to be below 0.01%, and the relative phase shift is 
detected within ±1.5°. No relative sign flips between the 
first- and second-harmonic amplitudes were detected. 



probe reg. 



C 



real 



Under present experimental conditions, this coefficient, 
Co, is found to be 1.028(3), and the measured PV param- 
eter is corrected accordingly. Principles of its derivation 
are given in Appendix [Cl 

In Fig. [Bl the PV interference parameter £//3 is shown 
as determined in 19 separate runs (^60 hours of integra- 
tion in total). Its mean value is 39(4) stat (3) syst . mV/cm, 
which is in agreement with the theoretical predictions 
[l5l |. The value of the PV parameter was extracted 
using the expression given in the first column of Ta- 
ble [Til taking into account the calibration correction 
C . Thus, |C| = 8.7 ± 1.4 x 10~ 10 ea , which is the 
largest APV amplitude observed so far (here we used 
|/3| = 2.24i°;?2 x 10" 8 ea /(V/cm) 0,0). 

The sign of the P V interference parameter £//3 is found 
by comparing the measurements with the theoretical 



VII. ERROR BUDGET 

The present measurement accuracy is not yet sufficient 
to observe the isotopic and hyperfine differences in the 
PV amplitude, which requires an accuracy better than 
w 1% for PV amplitude in a single transition fl7l - fl9l |. In 
the present apparatus the signal levels achieved values 
high enough to reach the signal-to-noise ratio (SNR) of 
2/^/t(s) for the PV asymmetry if the noise were domi- 
nated by the fluorescence-photon shot-noise (r is the in- 
tegration time). This is good enough to reach the sub- 
percent accuracy in a few hours of integration. However, 
a number of additional factors limit the accuracy. 

One of the most important noise sources is the fluctua- 
tions of the modulating- and DC-field parameters during 
the experiment. The first- and the second-harmonic sig- 
nals depend differently on the modulating electric field 
amplitude, Eq, and the DC-bias, thus, a noise in the 
electronics controlling the fields contaminates the first- 
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to-second harmonic ratio directly. A substantial effort 
was made to cope with this problem: a home-built high- 
voltage amplifier used in the first 13 runs was replaced 
by a commercial Trek 609B amplifier and a circuit con- 
trolling the DC-bias was upgraded. This allowed us to 
control the DC-bias and Eq with mV-scale accuracy that 
would make the SNR to approach the shot-noise limit if 
this were the only source of the noise. As seen in Fig. [5J 
the last six measurements exhibit higher accuracy than 
the rest. These are the runs after the HV-syst em u p- 
grade. However, the present SNR of s» 0.03/a/t(s) is 
worse by almost two orders of magnitude than the shot- 
noise limit. 

There are other fluctuations in the system parameters, 
such as light intensity fluctuations in the PBC, fluctua- 
tions of the spectral position of the PBC resonance with 
respect to the frequency reference, and noise in the de- 
tection system. All of them contribute to the noise in 
the first- and the second- harmonic signals but we found 
that such noise largely canceled in the ratio ru- 

However, there is a noise source, which is not canceled 
in the ratio. The following experiments demonstrated 
that this noise source is related to frequency excursions 
of the Fabry-Perot etalon serving as the frequency ref- 
erence for the optical system. In these experiments the 
excitation light was frequency tuned to a wing of the 
atomic resonance, and the first- and second-harmonic sig- 
nals were recorded without scanning over the resonance. 
Then, the same was done when the spectral position was 
set at the peak of the resonance, and a change of the SNR 
for the harmonics ratio was determined. These experi- 
ments were performed using the upgraded HV-system. 
Results of the measurements are presented in Fig. O For 
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FIG. 9: (color online) Impact of the frequency excursions of 
the Fabry-Perot etalon on the noise level in the harmonics 
ratio. A change in the noise level when the optical system 
was tuned from the wing of the atomic resonance to its peak 
is shown. In the inserts above the excitation light spectral 
position is shown schematically with respect to the atomic 
resonance. Arrows denote the fluctuations. 

a shot-noise-limied signal, the SNR at the peak of the 



resonance is expected to be a factor of about v2 higher 
than at the wing due to larger signal. It was found, 
however, that the SNR went up by a factor of 4 by tun- 
ing from the wing to the peak of the resonance. This 
demonstrates that the main source of noise is not photon 
statistics but fluctuations in the spectral reference. In- 
deed, the frequency excursions at the wing of the spectral 
line produces substantially more intensity noise due to a 
steeper spectral slope than that at the peak, where the 
slope is nominally zero. It must be emphasized, that in 
the case of slow frequency excursions (compared to the 
E- field modulation period), the noise in the first- and the 
second-harmonic channels would be canceled in the ratio. 
However, fast excursions can generate noise in the signal 
ratio. 

The factors affecting the measurement accuracy men- 
tioned above have an impact on the statistical error of 
the result. The present systematic errors (summarized in 
Table \TV\ has nearly the same significance as the statis- 
tical one and also comprises a number of factors. 

One of the most significant factors is the uncertainty 
in the field-imperfection contributions, Eq. (|2ip . This 
uncertainty is mostly due to statistical factors such as 
laser drifts, nevertheless, it provides an offset to the PV- 
parameter. Since the measurement of this contribution 
is actually the same measurement as the PV-effect, any 
improvements of the stability reduces the overall system- 
atic uncertainty. We would like to emphasize also that 
Eq. (|21[) represents a mean value of the imperfection con- 
tribution over numerous experiments averaging over pos- 
sible fluctuations of the field-imperfection contribution. 
These fluctuations may be partially responsible for the 
variance in the PV-parameter, and, thus, the statisti- 
cal uncertainty of its value. This fact demonstrates that 
the elimination of the field-imperfections is an essential 
requirement for improving the overall accuracy of the ex- 
periments. 

Another significant source of the systematic uncer- 
tainty is the uncertainty in the value of the electric field 
in the interaction region. While the voltage applied to 
the E-field plates and the correction electrodes is con- 
trolled precisely, the actual E-field value used in the PV 
parameter determination depends on the accuracy of the 
numerical modeling of the electric-field distribution in the 
particular geometry. There are two factors in the model 
contributing to the uncertainty: finite accuracy of mea- 
surements of the interaction region geometrical parame- 
ters, and the imperfect approximation of the geometry in 
the numerical simulation. 

However, while this systematic uncertainty plays a sig- 
nificant role for measurements of the PV parameter of a 
single isotope, for the isotope ratios this uncertainty will 
cancel (or will be substantially reduced), if the measure- 
ments observing different isotopes are performed without 
changing the E-field geometry. The same is true for the 
calibration parameter, Co, which also cancels in the iso- 
tope ratios. 

There are other, rather minor, factors contributing to 
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the systematic uncertainty, for example, a finite accuracy 
of the polarization angle flip, errors in the lock-in ampli- 
fiers, a finite dynamic range of the lock-ins etc. The net 
contribution of these factors is found to be < 1%. 

A summary of the systematic error budget is presented 
in Table |IV] 



TABLE IV: List of factors contributing to the systematic un- 
certainty of the PV parameter, 



Factor 


Uncertainty (%) 


E value: 




geometry 


5 


numerical modeling 


3 


E-field imperfections 


5 


Phase mixing 


0.5 


Other 


1 


Total (in quadrature) 


8 



VIII. TOWARDS MEASUREMENT OF THE 
ISOTOPE RATIOS 

As pointed out above, a goal of the future measure- 
ments of the parity-violation effects in ytterbium is ob- 
serving a difference in the PV effect between different 
isotopes. The net uncertainty of the PV parameter of a 
single isotope must be better than 1% based on the the- 
oretical predictions. To this end, a program of the appa- 
ratus upgrades and improvements is developed. Besides 
general improvements of the stability of the system pa- 
rameters, increase of the signal levels, suppression of the 
electronics noise etc., the main focus is on elimination 
of the frequency excursions of the frequency reference, 
which is a major source of the statistical noise. Improving 
the statistical uncertainty will contribute to more precise 
measurement and control of the E-ficld-imperfection con- 
tribution to the systematic part of the uncertainty. The 
latter is another high-priority improvement essential for 
reaching the goal. 

In the future apparatus, the referencing of the opti- 
cal system to the Fabry-Perot etalon will be replaced 
by locking the system to a femtosecond frequency comb 
that will be available shortly. The impact of the E-field 
imperfection is planned to be substantially suppressed 
by redesigning of the interaction region to provide more 
uniform and controlled electric field distribution. Until 
now, no scientific or technical obstacles were discovered 
preventing us to improve the apparatus to the desired 
level of sensitivity. 
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Appendix A: Derivation of transition amplitudes 

The total Hamiltonian (before including light-atom in- 
teractions and assuming B is along z ) can be written as 

H = H Atomic + -freeman + -^Stark + H^py, (Al) 

where H Atomic is the atomic Hamiltonian and -Hzeeman> 
H stark, and Hapv represent the contributions from the 
static magnetic field B, the static electric field E, and 
the parity non-conserving weak interaction, respectively. 
Here 



Zee man 



= —fi ■ B = .g/i B J • B = g/j, B J z B, 



(A2) 



where \x = —g^3 is the magnetic dipole moment of the 
atom, g is the Lande factor, /ib is the Bohr magneton, 
and J is the angular-momentum operator. Similarly, 

-ffstark = -d E = -diEi, (A3) 

where d is the atomic electric-dipole operator. Finally, 



H 



APV 



(A4) 



where Hq is a scalar operator. Summation over repeated 
indices is assumed. 

In the presence of a strong magnetic field, that is, 
when Zeeman splitting dominates Stark-shifts, it is use- 
ful to think of H\ = -ffstark + H^pv as a perturbation to 



H = H 
states 



Ator 
I2S+1 



nic + -Hzccman- In this case, the LS'-coupled 
Lj;M), such as | 3 D X ;M) and ^SojO), are 



eigenstates of the unperturbed Hamiltonian Hq- Then 
the first-order perturbation theory can be used to deter- 
mine the eigenstates of the total Hamiltonian: 



m =\a) + 2, -T~\ TTv 



(A5) 



where w(a) is the energy of state \a). (Perturbed eigen- 
states are denoted using an overbar.) 

The electric-dipole amplitude for the optical transition 
of interest is 



.4 



(El) 



m -( 3 D i; M|(-d.£)|iS )^ M 



(Stark) -(APV) 



where 



(A6) 



.4 



(Stark) _ 
M - 2 4 



E 



( 3 D 1 ;M|d-E| a ')(a'|d-,g| 1 S ) 
w( 3 Di) -u(a') 

( 3 n 1 ;M\d-e\a'}{a'\d--EfSo) 



(A7) 
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and 



Appendix B: Characterization of the PBC mirrors 



.4 



(APV) 
M 



E 



^DiiMliffgKX^Id-gpSo 
w( 3 Di) - w(a') 



£-f uj^So) - L>(a>) 

a' 

The Stark amplitude can be written as 

where = EiSj and 

T T _ di\a')(a'\dj dj\a'){a'\di 
13 ~ £-f cj( 3 D x ) - w(o') w( x S ) - w(o') ' 



(A8) 



(A9) 



(A10) 



Let and U q represent the irreducible spherical com- 
ponents of the tensors and EZy. Then TijUij = 
(-lyT^U* and Eq. jHI becomes 



A (st ai -k) = (_i)g T fc g< 3 Di . M |^*|i So) 



= (-1) 



k (3 Dl ||^||iS ) 



(0,0;fc,g|l,M) 



q 

= z/3(-l) 9 (E x £)_ 9 (0,0; 1,?|1,M). (All) 

Here /3 is the vector Stark transition polarizability and 
defined by 



1 



(A12) 



To derive Eq. (|ATT|) . we used (0, 0; fc, q\l, M) = S kl S qM 
zrxdTl q = J2 quq2 {l, qi ;l,q 2 \l,-q)E qi £ q2 = (i/^)(Ex 

For the parity-violating contribution to the El transi- 
tion amplitude we can likwise write 



4f V) = i£ i { 3 D 1 ;M\W i \ 1 S 



(A13) 



where 



H$\a'){a'\d, 



di\a')(a'\H° 



f w( 3 Di) - w(a') V w ( 1S o) - w(<»0 



A " ( A14 ) 



The finesse of the cavity is measured using the cavity- 
ring-down method • The laser beam is sent through a 
Pockels cell (Cleveland Crystals Inc. QX 1020 Q-Switch) 
and a polarizer before entering the cavity. The polarizer 
is aligned with the laser polarization so that the light is 
transmitted when there is no voltage applied to the Pock- 
els cell. A high-voltage pulse generator is used to send 
a fast step signal (30-ns wavefront) to the Pockels cell 
which rotates the polarization of the light so that it is not 
transmitted through the polarizer. The laser frequency is 
locked to the resonance frequency of the cavity, and then 
the Pockels cell is switched into the non-transmitting 
state, causing a fast interruption of the laser power. The 
subsequent decay of the light inside the cavity is moni- 
tored with a fast photodiode (50-MHz bandwidth) mea- 
suring the power transmitted through the back mirror of 
the cavity. The signal is fit to an exponential decay. The 
decay time is related to the finesse of the cavity (J 7 ) by 



T 



L 



T, 



where c is the speed of light, L is the cavity length, and 
r is the intensity decay time. An example of the PBC 
transmission signal and its fit are shown in Fig. 1101 
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Let £ q and W q represent the spherical components of 
the vectors £i and Wi, respectively. Then £iWi = 
{-l) q £i q W% and we have 



A% pv) =i(-iy£ 1 J 3 -D 1 -,M\W'\ 1 S ) 



= ^-l)^ g (3Dl|l ^ 1||1So) (0,0;l ig il,M) 



= »C(-l)«£l,<0 J 0;l 1 ?|l,M). 
Here £ is given by 



(A15) 



(A16) 



FIG. 10: (color online) Application of the cavity-ring-down 
method for the determination of the finesse of PBC. 

Following the analysis discussed in [2l[, if we denote 
the transmission of mirrors 1 and 2 by 1\ and T 2 , re- 
spectively and the absorption+scatter loss per mirror 
as h t 2 = (A + S)x t 2, then the total cavity losses C — 
Ti+T2+h + h determine the finesse T: 



T 



2n 



T 1 +T 2 + l 1 + l 2 



(Bl) 



Information on the transmission of the mirrors discrim- 
inated from the A+S losses can be obtained using the 
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measured value of the finesse and the power transmitted 
trough PBC, P tr : 



^ = 4TiT 2 ( — 



(B2) 



where Pi n is the input power, and e is a mode-matching 
factor. For two arbitrary mirrors, for which neither T 12 
nor lx t 2 are known the Eq. (|B1IB2I) do not provide a 
solution, since a number of variables exceeds the number 
of equations. Nevertheless, for two mirrors from the same 
coating run when one can assume that 1\ = T 2 = T and 
h = h = I, the equations ()B1IB2[) become 



T = 



T + l 



Ptr 



AT 2 



and for known mode-matching factor e the parameters of 
the mirrors (T and I) can be determined. The factor e 
depends on the geometry of the cavity, and is assumed 
to stay constant upon replacing of the mirrors, if the ge- 
ometry of the input laser beam and the configuration of 
the PBC are unchanged. This gives the possibility to 
calibrate this factor by using a mirror set for which the 
transmission is known. We used for this purpose the mir- 
ror set purchased from Advanced Thin Films, Inc., for 
which reliable data on the transmission of the mirrors 
is provided by the supplier. By measuring the finesse 
of the PBC comprised of these mirrors and the ratio of 
the transmitted-to-input power, the mode-matching fac- 
tor and the A+S mirror losses I are found. This set is 
not an actual mirror set that was used in the PV experi- 
ment, nevertheless, the parameters of other mirrors were 
determined by replacing one mirror in the "reference" 
set by the "test" mirror, parameters of which are sought. 
The geometry of the cavity was unchanged during the 
replacement. This tactic allows for the measurement of 
parameters of any arbitrary mirror. 



Appendix C: Impact of the phase mixing effect on 
the harmonics ratio 



Atoms undergo the 6s 2 1 Sq — > 5d6s 3 D! transition 
in the interaction region where they are illuminated by 
408- nm light and are exposed to the static magnetic field 
and the oscillating electric field E(t). Excited atoms 
then spontaneously decay from the 5d6s 3 D! state to the 
metastable 6s6p 3 Po state. The population of 6s6p 3 Po 
is proportional to the transition rate Rm for M = 0, ±1. 
Without loss of generality, we assume that the constant 
of proportionality is equal to one. 

The rate Rm is measured in the probe region. The 
probe region is located a distance d sa 20 cm away from 
the interaction region. Therefore, an atom that arrives 
at the detection region at time t experienced an electric 
field with magnitude E(t—d/v z ) in the interaction region, 



where v z is the atom's speed and d/v z is the amount of 
time required for the atom to travel a distance d. 

Because some atoms travel faster or slower than oth- 
ers, the detection region is full of atoms that have each 
experienced a different electric field while in the inter- 
action region. Each atom contributes to the total rate 
and hence the observed rate Rm is the thermal average 
of every contribution: 



R M (t;ui,d,v ) = / R M (t - d/v z )f(v z ;v )dv z , (CI) 
Jo 



where 

/(«,; wo) dv z = 2(v z /v ) 3 e-^/ v ^dv z /v , (C2) 

is the probability for an at om to hav e speed between v z 
and v z + dv z . Here v = ^/2~k^f~fm = 2.9 x 10 4 cm/s is 
the thermal speed, T w 873 K is the oven temperature, 
and m = 161 GeV/c 2 is the atomic mass of Yb. 

It is convenient to introduce the dimensionless vari- 
ables x — v z /vq and t = Lit, and the dimension- 
less parameter a = uod/vQ. Then the average rate 
RM{t;uj,d,vo) — > Rm{t;oi) depends only on the dimen- 
sionless quantities a and r, and Eq. (|C1[) becomes 



R m (t- a) = 1ZZ + n M \I(a)\cos( T + Arg[7(a) 

,[2] 

with 



1Z l *}\I(2a)\ cos(2r + Arg[I(2a)]), (C3) 



1(a) 



e~ la/x f{x-l)dx. 



(C4) 



Note that \I(a)\ — > as a — > oo whereas \I(a)\ ~ 1 when 
a < 1. This places a limit on the modulation frequency: 
We require that oj < vo/d = 2ttx 230 Hz in order to 
avoid a significant decrease in signal. 

The lock-in amplifier receives an input signal propor- 
tional to Rm and returns two output signals and S^j 
corresponding to the first and second harmonic compo- 
nents, respectively. This process can be modeled as 



S 



M 



i>n, a) = 



1 



Rm (t; a) cos(nr + <fi n ) dr 



TZ M ] \I{na)\ cos(Arg[/(na)] + </)„), (C5) 



where the phases ipi^ of the lock-in amplifier are chosen 
to maximize the signals S M ' . That is, 

4> n = <j) n {a) = -Arg[/(na)]. (C6) 

Our measurement sm is the ratio of the first- and second- 
harmonic signals: 



sm 



S^faa) TZ [ lj\I(a)\ 



S M l (<f> 2 ;a) KW\I(2a 



r M xC(a), (C7) 
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where C(a) = |/(a)|/|/(2a)| is the correction factor. 
Therefore, we must further divide the ratio sm of ob- 
served output signals by C(a) to measure the ratio r M . 

The correction factor C(a) and the optimal lock- in 
phases ^1,2(0:) inherit dependence on the modulation fre- 
quency (u> = 2tt x 76.2 Hz), the distance between inter- 
action and detection regions (d ft 20 cm), and the oven 
temperature (T ft 873 K) through the parameter a: 



>d 



a 



0.33(2), 



^2k B T/m 
where the uncertainty in a is given by 



da = ay / (ST/2T) 2 + (5d/d) 2 , 



(C8) 



(C9) 



for ST « 50 K and 5d ft 1 cm. The correction factor can 
be computed numerically and has a value 



C = C(a) = 1.028(3), 



(CIO) 



with uncertainty given by 5Co = \C'(a)\Sa. Likewise, 
the lock-in phases have the following values 



= Ma) = 16(1)°, M> = <hi*) = 33(2)°, (Cll) 



where 5<fi n o = \<f)' n (a)\ 5a. 

In order to understand the impact of imperfect phase 
selections, we include the effects of slight deviations from 
the optimal phase <j> n {oi) by taking 



K -> 4>n{a) + Ifn, 



(C12) 



where ip n ft represents a small deviation. Then the 
correction factor becomes 



C(a) -> C(a, <p u ip2) = C{a) x c ^pl ( C 13) 



C0S(^2 



and hence Co = C(a, 0,0) = C(cv) = Co- The uncer- 
tainty in the correction factor becomes 



<5C -> (5C - ^8C% + 8<p\ + Si4, (C14) 



where <5<^„ is the uncertainty in the deviation (£>„. To de- 
rive this expression, we estimated the partial uncertainty 
in C due to ip n by d* n C{a, <p\,yi) Sip 2 n . 

To estimate the uncertainty 6<p n , we assume that we 
are within about 1° of the optimal phase. This choice 
is consistent with the magnitude of the uncertainty in 
the optimal phases <pio and (f>2o- Therefore, we will take 
5ip n = 5<p n o to be the accuracy with which we can select 
the lock-in phases. Then 5ipi = 0.02, 5(p 2 = 0.03, and 



5C = 0.0031 « 0.0029 = SC . 



(C15) 



Hence small deviations (on the order of 1°) have a negli- 
gible effect on the uncertainty in the correction factor. 
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